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Abstract 

We study the evolution of parton distributions at small x values, obtained in 
the case of flat initial conditions. The contributions of twist-two and (renormalon- 
type) higher-twist operators of the Wilson operator product expansion are taken 
into account. The results are in good agreement with deep inelastic scattering 
experimental data from HERA. 



1 Introduction 

The measurements of the deep- inelastic scattering structure function (SF) F2 in HERA 
m 12] have permitted the access to a very interesting kinematical range for testing the 
theoretical ideas on the behavior of quarks and gluons carrying a very low fraction of 
momentum of the proton, the so-called small x region. In this limit one expects that non- 
perturbative effects may give essential contributions. However, the reasonable agreement 
between HERA data and the next-to-leading order (NLO) approximation of perturbative 
QCD that has been observed for > IGeV^ (see the review in Ref. indicates that 
perturbative QCD could describe the SF evolution up to very low values, traditionally 
explained by soft processes ^. It is of fundamental importance to find out the kinematical 
region where the well-established perturbative QCD formalism can be safely applied at 
small X. 

The standard program to study the small x behavior of quarks and gluons is carried 
out by comparison of data with the numerical solution of the DGLAP equations by fitting 
the parameters of the x profile of partons at some initial and the QCD energy scale A 
(see, for example, [HI El)- However, if one is interested in analyzing exclusively the small 
X region {x < 0.01), there is the alternative of doing a simpler analysis by using some of 
the existing analytical solutions of DGLAP in the small x limit (see, for example, Ref. |Sj 
for review). This was done so in Refs. [H1-J2] where it was pointed out that the HERA 
small X data can be interpreted in terms of the so called doubled asymptotic scaling 
phenomenon related to the asymptotic behavior of the DGLAP evolution discovered in 
[THj many years ago. 

Here we compile results obtained in 0-^1], where the contributions of higher- twist 
(HT) operators (i.e. twist-four ones and twist-six ones) of the Wilson operator product 

^The agreement can be explained 0j by an effective coupling constant scale essentially higher then 
(see review and references therein). 



expansion are taken into account ^. The importance of the contributions of HT operators 
at small-x has been demonstrated in many studies (see fT3]). 

We would hke to note that the results of are the extension to the NLO QCD 
approximation of previous leading order (LO) studies [THl 15]. The main ingredients are: 

1. Both, the gluon and quark singlet densities are presented in terms of two compo- 
nents ('+' and '— ') which are obtained from the analytical dependent expressions of 
the corresponding ('+' and '— ') parton distributions moments. 

2. The '— ' component is constant at small x, whereas the '+' component grows at 
> Qo as ~ exp [o^lo)-, where 

(yNLO = 2\/ {d+s + D+p) Inx, 

and the LO term d^ = — 12//3o and the NLO one = + d^Pi/Po with 6/++ = 
412//(27/?o)- Here the coupling constant = as/(47r), s = \n{as{Ql)/as{Q^)) and 
P = CLsiQo) ~ (^siQ"^), Po and Pi are the first two coefficients of QCD /^-function and 
/ is the number of active flavors. 



2 Basical formulae 

Thus, our purpose is to demonstrate the small x asymptotic form of parton distributions 
in the framework of the DGLAP equation starting at some Qq with the flat function: 

/a^(Qo) = ( hereafter a = (1) 

where /^^ are the leading-twist (LT) parts of parton (quark and gluon) distributions 
multiplied by x and Aa are unknown parameters that have to be determined from data. 
Through this work at small x we neglect the non-singlet quark component. 

We would like to note that new HERA data [2] demonstrate a rise of SF F2 at low 
values {Q^ < IGeV^) when x ^ (see Fig. 2, for example). The rise can be explained in 
natural way by incorporation of HT terms in our analysis (see section 2.2). 

We shortly compile below the main results found in ^ ^] at the LO approxima- 
tion (the LT results at the NLO approximation may be found in [0]). The full small x 
asymptotic results for parton distributions and SF F2 at LO of perturbation theory is: 

F^ix, Q') = e- [/f (x, Q') + (x, Q') + (x, Q')] , (2) 

where e = (J^i^D/f is the average charge square. The LT parts /^^(x, Q^) and the HT 
ones fl^'^{x, Q"^) can be represented as sums of the '+' and '— ' components 

fti^,Q') = lT'H^,Q') + fT'-i^,Q')- (3) 



2.1 The LT contribution 

The small x asymptotic results for parton distributions have the form 

/f'+(x,g^) = (A, + ^A,)/o(a)e-^+«^ + 0(p) , (4) 

•^Refs. fTUllll) and the present note deal with the contributions of HT operators only in the renormalon 
model approach (see JJ). The complete analysis of the HT contributions to the SF F2, dF2/d\n{Q'^) 
and will be presented in |12|. 
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n'^H^,Q') = ^fA, + ^AJp/i(a)e-'^+«^ + 0(p), (5) 



9V " 9 

/f'-(x,g2) = ~A,e-''~^'> + Oix), (6) 
/f'-(x,Q2) = A,e-'-^'> + 0(x), (7) 

where = 1 + 20//(27/?o) and = 16//(27/3o) are the regular parts of d+ and 

(i_ anomalous dimensions, respectively, in the limit n ^ 1 ^. The functions I^, {v = 0, 1) 
are related to the modified Bessel function I^, and to the Bessel function by: 



The variables a and p are given by 



p'luicr), if s > 

{—'pY-Jyia), p = ip, a = ia, if s < 



a = 2^d,sHx), P=^^^ = ^i^ (9) 

2.2 The HT contribution 

Using the results in [IHIIII] (which are based on calculations IHIIIH]), we demonstrate 
the effect of HT corrections in the renormalon model case (see recent review of renormalon 
models in jT7|). We sketch below the basical results making the following substitutions 
in the corresponding twist-two results presented in Eqs.(jlI)-(IZj) (complete formulae can be 
found in |I2!): 

/;2.+ (a;,Q2) and f^^'+{x,Q^) (see Eqs.dH) and ®) ^ f!;^'+{x,Q^) and f^^'+{x,Q^) by 

AM.) ^ ^„ W 1^ ,„.(^ _ ,„ . ,10) 

where fci = 1, ^2 = —8/7 and Ai^^ and A2,a are magnitudes of twist-four and twist-six 
corrections. 

/f.-(x,g2) and '-(x,Q2) (gee Eqs.dHD and ©) ^ f;^'-{x,Q^) and f^^'-{x,Q^) by 

From Eqs. (jl(jp - (jl3|) one can notice that the higher- twist terms modify the flat condition 
Eq.(IT)). They lead to a rise of parton distributions and, thus, SF F2 at low value Ql, when 
X — > 0. This is in excellent agreement with the recent low HERA data [2]. 



■^From now on, for a quantity k{n) we use the notation k(n) for the singular part when n ^ 1 and 
k{n) for the corresponding regular part. 
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Fig. 1 





Q'=1Z GeV 





Q"' = 6.5 GeV 






" = 8.5 GeV 




Q'' = ZO GeV 




X 



X 



Figure 1: The SF F2 as a function of x for different bins. The experimental points are 
from HI pp. The inner error bars are statistic while the outer bars represent statistic and 
systematic errors added in quadrature. The dashed and dot-dashed curves are obtained 
from fits (based on LT formulae) at LO and NLO respectively with fixed Qq = 1 GeV^. 
The solid line is from the fit at NLO giving Ql = 0.55 GeV^. 



3 Results of the fits 

With the help of the results obtained in the previous section we have analyzed F2 HERA 
data at small x from the HI and ZEUS collaborations Plf2]. In order to keep the analysis as 
simple as possible we have fixed the number of active flavors /=4 and ^j^|g"(^/ = 4) = 250 
MeV, which is a reasonable value extracted from the traditional (higher x) experiments. 
Moreover, we put Ai^^ = A2,a in agreement with ^H]- The initial scale of the parton 
densities was also fixed into the fits to Qo = 1 GeV^, although later it was released to 
study the sensitivity of the fit to the variation of this parameter. The analyzed data region 
was restricted to x < 0.01 to remain within the kinematical range where our results are 
accurate. 

Fig. 1 shows F2 calculated from the fit (based only on LT formulae) with > 1 GeV^ 
in comparison with 1994 HI data (first article in P). Only the lower bins are shown. 
One can observe that the NLO result (dot-dashed line) lies closer to the data than the LO 
curve (dashed line). The lack of agreement between data and lines observed at the lowest 
X and bins suggests that the fiat behavior should occur at lower than 1 GeV^. In 
order to study this point we have done the analysis considering Qq as a free parameter. 
Comparing the results of the fits (see [H]) one can notice the better agreement with the 
experiment at fitted Qq = 0.55 GeV^ (solid curve) is apparent at the lowest kinematical 
bins. 

Fig. 2 shows F2 calculated from the fit at LO (based on LT & HT formulae) in 
comparison with 1995 HI and ZEUS data |2|. One can observe that these results (solid 
line) lies closer to the data than the twist-two results (dashed line). We have done the 
analysis considering Qq as a free parameter. Comparing the results of the fits (see [T^ ) 
one can notice the better agreement with the experiment at fitted Qq = 0.61 GeV^, which 
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Fits to H195 data. Also shown ZEU95 data. dashes = LO solid = LO + HTren. nf = 4 
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Figure 2: The SF F2 as a function of x for different bins. The experimental points are 
from HI and ZEUS [2]. The inner error bars are statistic while the outer bars represent 
statistic and systematic errors added in quadrature. The solid curves are obtained from 
fits at LO, when contributions of HT terms have been incorporated. The dashed curves 
show only twist-two contributions. 

is close to Ql in the analysis of 1994 HI data (see Fig. 1). 



4 Conclusions 

We have shown that the results developed recently in 0- [12] have relatively quite simple 
form and reproduce many properties of parton distributions at small x, that have been 
known from global fits. 

We found very good agreement between our approach based on QCD and HERA 
data, as it has been observed earlier with other approaches (see the review [3_.). The 
(renormalon-type) higher-twist terms lead to the natural explanation of the rise of F2 at 
low X for the lowest values of (< 1 GeV^). The rise has been discovered in recent 
HERA experiments P|. 

As next step of our investigations, we plan to study contributions of HT operators to 
relations between parton distributions and deep inelastic structure functions, observed, 
for example, in [THl 120] • 
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